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Measurements of the rf penetration depth A(r, H, 9) are used to study the superconducting 
order parameter, vortex dynamics in the mixed state and dehneate critical fields in the boro- 
carbide superconductor YNi2B2C. The lower critical field has an anomalous T dependence, 
Hci{T) — 1.12 [1 — (T/Tc)] kOe, which is however consistent with independent superfluid den- 
sity measurements at microwave frequencies. The vortex response is dominated by viscous flux 
flow, indicative of extremely weak pinning, and is parametrized by a field scale Hc2.cS- The angu- 
lar dependence of the vortex contribution \(6) is in good agreement with the Coffey-Clem model. 
Structure is seen in the depairing transition in the vicinity of the upper critical field, with the exis- 
tence of well-defined critical fields Hc2a, Hc2b and Hc2c, with the vortex field scale Hc2,cS closest to 
Hc2b- Overall the measurements indicate that YNi2B2C has a rich and unusual field dependence of 
its transport parameters. 
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The borocarbide superconductors LnNi2B2C where Ln 
= {Y, Lu, Tm, Er, Ho and Dy} [0-^ are a recently 
discovered family of superconductors with elevated tran- 
sition temperatures. While there have been numerous 
studies of these superconductors using traditional probes 
such as magnetization, specific heat, tunneling, etc. there 
have been relatively few which directly probe the order 
parameter and dynamics of vortices 

In this paper we report on studies of the radio fre- 
quency (rf) penetration depth X{T, H, 0) of single crys- 
tal YNi2B2C (Tc — 15.5 K) which were carried out in 
a highly sensitive rf tunnel diode oscillator set-up. The 
sample is placed in an inductive coil which forms part of 
a tank circuit, which is self-resonant typically at 3 MHz 
and is driven by a tunnel diode circuit. Changes in the 
penetration depth or skin depth of the sample, caused by 
varying temperature, T , applied dc magnetic field, H , or 
angle between H and the ah plane, 6, lead to changes in 
the coil inductance. These are then detected as changes 
in resonant frequency, and converted to changes in A 
via the relation 5\{T, H, 0) = -g [f{T, H, 9) - /o(T, H)], 
where g is a geometrical factor set by the sample dimen- 
sions and /o refers to the resonant frequency without the 
sample. The very high stability of the circuit, typically 
1 Hz in 3 MHz, leads to a very high resolution of the 
order of a few A. The technique has been extensively 
validated through precise measurements in the cuprate 
superconductors of the non-linear Meissner effect and of 
vortex parameters such as Hd and pinning force con- 
stants In one experimental set-up the sample can 
be oriented with H \\ c {9 = 90°) or H ± c {9 = 0°), and 
fields up to 7 T can be applied. In another set-up the 



angle 9 can be varied continuously between 0° and 90° 
with an angular accuracy of 0.2°. Here the maximum 
applied field is 6.4 kOe. In all cases the rf field Hii \\ ah 
and i?rf -L H . 

The single crystals of YNi2B2C used for these mea- 
surements were grown using the Ames Lab Ni2B fiux 
growth method |^,|^. Crystals grown via this method 
are in the form of plates with the crystallographic c axis 
perpendicular to the surface of the plate. The crystal 
used for these measurements has approximate dimensions 
1.4 X 1.2 X 0.2 mm^. 



I. LOWER CRITICAL FIELD: Hci 

Typical results for A\{T,H,9) = 5\{T,H,9) - 
6X{T, iJ = 0) vs. H are shown in Fig. |, for both 6* = 0° 
and 90°. The low field portion is further shown in de- 
tail in the inset, which reveals essentially no change until 
a critical field Hd is reached, above which A increases 
rapidly. This critical field Hd represents the field at 
which fiux first enters the sample and is governed by 
both the Meissner state of the superconductor and an 
effective geometric barrier at the surface. We are justi- 
fied in calling it the lower critical field as demagnetiza- 
tion corrections are negligible in the case where H \\ ah. 
Note that the presence of surface barriers would lead to 
hysteresis in the field dependence of AX{T, H,9) at the 
onset of fiux entry. In our data, the observed hystere- 
sis is minimal and is within the experimental resolution. 
It is also to be noted that surface barrier effects would 
scale the field values but do not affect the temperature 



1 



dependence. 




low Tc superconductors and even differs from the behav- 
ior of YBa2Cu307_5. However this anomalous behavior 
finds support in measurements on a related quantity, the 
superfluid density, measured via the condensate contri- 
bution to the complex conductivity, a2, which we have 
measured elsewhere Q using a 10 GHz Nb cavity. There, 
(7i was shown to have an anomalous temperature depen- 
dence well represented by o^iT) = 10^(1 — t) (Om)~^ 
(also shown in Fig. ||). Now (T2 ~ I/^iqIjjX^, and since 
Hci ~ ($o/47rA^) [ln(K) -1-0.5], there is a good correla- 
tion between the linear temperature dependence of the 
Hci and (T2 measurements. (Note that for clean type-H 
superconductors, k varies as 1/(1 -\-t^) and thus ln(K) is 
expected to have a weak T dependence 0). 



II. VORTEX DYNAMICS AND DEPAIRING 



FIG. 1. AA ?;s. // at T = 8 K for YNiaBaC. (Inset) Ex- 
panded plot of low field portion, showing onset at Hd- A 
schematic of the field/sample orientation is also presented as 
an inset 

Notice that unlike magnetization measurements, where 
Hci is deduced from deviations from linearity in the 
M{H) curve, the signature in our data at Hd is very 
sharp. This is because the experiment effectively mea- 
sures the flux density B{H) (as will become evident 
later). Hence, the signature at Hd is a change from 
zero, and not a change of slope as in M{H) data. 




FIG. 2. Hci vs. T. 02 vs. T from Ref. |g| is shown for 
comparison. The dotted line represents a BCS temperature 
dependence of Hd while the dashed line is a 1 — t dependence 
obtained from a least squares fit to the data. 

The T dependence of Hci{T) for YNisBaC is shown in 
Fig. ^. This essentially has a linear behaviour from Tc 
to 2 K, weU represented by HdiT) = 1.12(1 - t) kOe, 
where t = T /Tc. This is unlike that seen in conventional 



Above i?ci, X{T,H) increases strongly in the 9 = 0° 
case, due to the increasing B in the sample. A simi- 
lar increase also occurs in the 9 = 90° data, although 
here the effective Hd — > due to demagnetization. This 
field dependence X{H) for H > Hd arises from two ef- 
fects of the magnetic field on the superconductor: depair- 
ing of condensed electrons (which dominates when the rf 
current J^f || H ) and dynamic vortex response (which 
dominates when Jj-f _L -ff ). This is followed by a sharp 
increase starting and terminating at field values which we 
call Hc2a and Hc2b respectively. Above Hc2b, X increases 
almost linearly with H , which is consistent with the nor- 
mal state showing a strong positive magnetoresistance. 
There is an additional field Hc2c just above Hc2b which 
can be identified as the point where the 6* = 0° and 90° 
curves meet as seen in Fig. 1. 

We have observed a similarly rich structure in the tran- 
sition region of the Ln = {Er, Ho and Dy} compounds 
|llo| . The presence of multiple structure in the transition 
region, although interesting in its own right, makes the 
identification of a single upper critical field Hc2 some- 
what ambiguous. This difficulty is also compounded by 
the linear- _ff magnetoresistance in the normal state. All 
three field scales {Hc2a, Hc2b and Hc2c) agree quantita- 
tively with the Hc2 values reported in the literature 
While the exact origin of this structure is not yet precisely 
known, its presence in all the borocarbide superconduc- 
tors we have studied seems to suggest an intrinsic phe- 
nomenon and is not likely to be due to spurious chemical 
phases, as will be presented in a forthcoming publication 

In this paper, we focus on the general nature of the 
depairing and vortex contributions to the rf penetration 
depth. To test the validity of the standard Coffey-Clem 
model for the depairing and vortex response in this sys- 
tem, an effective upper critical field Hc2 coinciding with 
Hc2b is assumed for subsequent analysis. It is to be 
pointed out that our choice of Hc2b as the effective upper 
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critical field is not arbitrary and is based on other factors 
like the validity of the Bardeen-Stephen approximation 
PI which relates the upper critical field to the fiux flow 
resistivity. A discussion on fiux viscosity presented later 
in the paper will clarify this issue further. 



III. DETAILED COMPARISON TO THEORY 



We now turn to an understanding of the field depen- 
dence of X{T,H,9). The principal causes for this field 
dependence are the depairing effect of B{H) on the con- 
densate and the dynamic response of the vortices to the 
Lorentz force induced by the rf current. These contribu- 
tions have been treated self-consistently by Coffey and 
Clem, leading to the following expression for the com- 
plex penetration depth which includes the full angular 
dependence: 
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where Xl{T,B) and Snf{T,B) are the condensate pene- 
tration depth and the normal-fiuid skin depth in the su- 
perconducting state, and K,p{T) and r]{T) are the vortex 
pinning force constant and viscosity respectively. The 
field-dependence of the condensate background is rep- 
resented by Ai(r,B) = Aio/(T,B) and <5^/T,B) - 
Sl{T, B) [I - /(T, B)], where for instance /(T, B) = [l- 
[T/T^Y] [1 - {B/Bc2)\ in a two-fiuid model. For 6* = 0°, 
the rf currents in the ab plane are parallel to the vortices 
and hence the Lorentz force is zero for those vortices 
at the sample faces. (A schematic representation of the 
relative orientations of H^f and the rf current, Jrf, with 
respect to H is shown as an inset in Fig. ^ . Thus the vor- 
tex response is mostly absent in the 9 = Q° case, except 
for a small contribution of order 0{d/W) ~ 0.1 from the 
edges, where d and W are the sample thickness and width 
respectively. This can be incorporated phenomenologi- 



cally in terms of an aspect ratio r = 
which mixes the vortex contribution 
force-free {9 
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The angular 



0°) term as ^A^| r\\. 

dependence provides a nice way of separating out the de- 
pairing and the vortex contributions to the penetration 
depth. 




FIG. 3. 2'^^^A\{T,H)/S„{T,H) vs. H at T = 8 K and 
9 = 0° and 90°. The solid lines represent theoretical cal- 
culations based on Eq. |l] for the two orientations presented, 
including finite aspect ratio effects for the ^ = 0° case with 
r = 0.11. The dotted line ignores the effects of finite aspect 
ratio. The dashed horizontal line represents the normal state 
skin depth. 

In Fig. 1^ we present the same data as in Fig. |l| nor- 
mahzed with respect to 5„(T, _ff)/-\/2. Motivated by 
the observation that the field dependence of Sn{T < 
Tc,H > Hc2) follows that of (5„(T > T^, H) we have 
used Sn{T,H) « Sn(T > Tc,H) for this normalization. 
Theoretical calculations based on Eq. |l|, taking into ac- 
count the small aspect ratio as discussed earlier and also 
assuming negligible pinning (i.e. Kp — > 0), are shown 
in Fig. y along with the data. For the vortex contribu- 
tion, at fixed frequency / the response is pinning domi- 
nated if / < /c and viscosity dominated if / > /c, where 
fc = Kp/2TTri. It is quite striking that pinning is negli- 
gible even at 3 MHz, which would effectively imply that 
fiux is in free fiow as soon as vortices enter the sample 
above Hd- This is in contrast to the estimated fc values 
of 20 GHz in YBaaCuaOr-i ^ and - 100 MHz in low 
Tc materials jl^. For 9 = 0°, the low field variation of 
A is adequately reproduced by the calculation, confirm- 
ing the presence of a finite vortex contribution consistent 
with the aspect ratio independently determined from the 
sample dimensions. 

The 9 = 90° data shows a dip at a field below the effec- 
tive Hc2 ■ This can result from sudden changes in pinning 
or scattering which are parametrized by the terms Kp and 
1] in Eq. |l]. A rise in Kp at a threshold field can be asso- 
ciated with the "peak effect" commonly seen in type-II 
superconductors with weak pinning. The signature of the 
peak effect is a peak in critical current density and would 
show up as a dip in the rf susceptibility, which is what we 
are effectively measuring in the 9 = 90° case. There is a 
general consensus that the peak effect arises due to en- 
hanced pinning as a result of a rapidly decreasing rigidity 
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of the vortex lattice just below Hc2 flS)- However, there 
is still some controversy over whether it is a softening 
of the tilt modulus (C44) or shear modulus (Cee) that 
leads to this decrease in rigidity ||T^-[T^. An alternate 
possibility is that this dip feature is due to an anomalous 
field dependence of rj arising as a result of unusual quasi- 
particle scattering in the superconducting state. Since 
the effect is much stronger in the vortex response rather 
than the depairing, this is an indication that the scat- 
tering due to bound quasiparticle states in the normal 
vortex cores could also be responsible for the observed 
dip in the 9 = 90° data below the effective Hc2- 

When the vortex term dominates as at = 90°, 
the A(T, H) data can be used to extract the pinning 
forces and viscosity as we have done elsewhere in both 
single crystals and films of YBa2Cu307^5 ^ 
Bardeen-Stephen estimate of the viscosity can be made 
using VBsiT,B) = $oSc2(T)/p„(T,S), which yields 
??Bs(0,0) ~ 1 X 10-6 Ns/m^ for YNiaBaC, compared to 
values of 10"^ Ns/m^ for Nb. 




/i.H [T] 



the observed response and theoretical Bardeen-Stephen 
response (indicated by the solid line in Fig. ^) can be 
attributed to an underestimate of p„ (T, B) in the super- 
conducting state as well as some variation of Kp{T). If 
the dip observed in the 9 — 90° data is associated with 
a change in pinning, we can estimate the change in Hp in 
the peak region from the relative depth of the dip since 
VeS ~ ri/{l + Kp/ur]). A plot of /c,pcak = Kp/2TTr]{T) 
extracted in this form is shown in the inset to Fig. ^. 

IV. ANGULAR DEPENDENCE 
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FIG. 5. A\ vs. e at T = 8 K and H = 1 kOe. The line 
represents a least squares fit to Eq. (Inset) vs. sin^ 9 
for the same data for comparison with Eq. 0. A plot for poly- 
crystalline Nb at r = 7 K and H — 1.5 kOe is also shown. 



FIG. 4. Flux flow resistivity Pf/pn obtained from the 
9 = 90° data in Fig. ^ The solid line indicates the expected 
Bardeen-Stephen behavior. The dashed horizontal line indi- 
cates the normal state resistivity, (inset) fc = Hp/ri vs. T 
extracted from the relative peak height is presented as an 
inset. The horizontal line in the inset indicates the measure- 
ment frequency. 

The data in Fig. ^ can also be presented as the nor- 
mahzed flux flow resistivity p//p„ = 2X'^{T, B)/Sl{T, B) 
which is presented as a function of H for T/Tc = 0.5 
in Fig. |4| The flux flow resistivity approaches the nor- 
mal state almost linearly with a slope 1 / Bc2, cff (~ 1/ Bc2) 
consistent with the observations of Gittleman and Rosen- 
blum in conventional superconductors p^ . Overall the 
response appears to be well described by Eq. ^ assuming 
77 = riBsiT^B) and Hp — s- 0, consistent with magneti- 
zation experiments which have also shown evidence for 
extremely weak pinning |l9[. The differences between 



The full angular dependence of A for Hci < H < Hc2a 
is shown in Fig. || along with the fit generated using Eq. [|. 
The inset shows plotted against sin'^ 9 for the same 
data. A similar plot for a poly-crystalline Nb sample of 
comparable dimensions is also shown in the panel. The 
excellent agreement between the experimental data and 
the Coffey-Clem expression of Eq. |] is obvious from the 
good fit and also the linear dependence of A'^ with sin^ 9. 
The slopes are vastly different for YNi2B2C and Nb as 
the response is pinning dominated in Nb whereas pin- 
ning is almost irrelevant in YNi2B2C at 3 MHz. The 
Coffey-Clem calculation as presented in Eq. |l] assumes 
an isotropic Hc2 case and from the good fit with the 
YNi2B2C data one can infer that the effective Hc2 is gen- 
erally isotropic in this system Q , although a slight devi- 
ation of the data from the model at low angles might in- 
dicate the presence of some degree of anisotropy in these 
layered systems. Elsewhere we have observed unusual 
angular dependence and anisotropy effects in HoNi2B2C 
related to magnetic order and details about this will be 
given in a separate publication. 
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FIG. 6. Hc2a, Hc2b and //c2c '^s- T obtained from the 
6 — 0° data. Lines representing a (1 — depen- 
dence are overlaid. Hc2 values from Ref. Q are indicated by 
the thick dashed line. The region where the dip is observed 
in the 9 — 90° data is indicated by the dotted hatch lines. 

The magnitude and T dependence of the three features 
in the vicinity of iJc2 for YNi2B2C {Hc2a, H^jib and Hc2c) 
are shown in Fig. ^ are in good agreement with those 
determined by magnetization measurements |^ and the 
overall temperature dependence is fit reasonably well by 
a(l — form. Two aspects make the observed 

onset field i?c2a extremely interesting: Hc2a is located at 
about the same field where the dip is seen in the {9 — 90°) 
data and occurs further from the effective upper critical 
field than would otherwise be expected. The theoretical 
calculation fails to reproduce the finite slope of the tran- 
sition, which would imply an unusual field dependence of 
A(-B). Looking at it in a different way, the simple varia- 
tion of the superfluid density (= XfjX^) as (1 — B/Bc2) 
within the framework of a two-fluid model breaks down 
and one has to assume a stronger reduction in the su- 
perfluid density as B Bc2- Note that earlier studies 
on YNi2B2C 1^ using a superconducting Nb cavity at 
10 GHz also indicated disagreement with BCS or two- 
fluid models in the T dependence of the surface resistance 
and penetration depth, as displayed in Fig. ^ in terms of 
tT2(r). This may indicate the existence of scattering of 
unknown origin in the superconducting state. We be- 
lieve the origin of Hc2a could result from this anomalous 
scattering or vortex lattice softening/melting or a com- 
bination of both these effects. Whatever the cause, the 
underlying dependence of AA(T, H,6 = 0) in the vicinity 
of Hc2a is masked to some extent by the mixing of vortex 
and depairing terms in the (9 = 0°) data. 

In conclusion, the rf measurements presented here re- 
veal a rich variety of electrodynamic response in the su- 
perconducting state of the non-magnetic borocarbide su- 
perconductor YNi2B2C. In particular, \{T,H = 0) and 
Hci both have anomalous non BCS T dependences. The 



experiments have very high sensitivity which enable ob- 
servation of subtle effects on the superconducting order 
parameter and the vortex lattice due to the applied field. 
A unique way of studying the vortex response and field- 
induced depairing individually by variation of 9 is pre- 
sented. The angular dependence of A(T, H, 9) is well de- 
scribed by the Coffey-Clem theory with an isotropic Hc2 ■ 
A dip is observed in the vortex response which can be 
related to a change in pinning (and thus to the "peak 
effect") or alternatively to anomalous scattering in the 
mixed state. Pinning is found to be extremely weak, and 
consequently, the crossover frequency is expected to be 
orders of magnitude lower than that of most conventional 
and high- Tc superconductors. Structure is observed in 
the depairing transition from the superconducting to nor- 
mal state, which is also found to have an anomalous H 
dependence. Thus although YNi2B2C is in many ways 
a conventional high-K type-II superconductor, the high 
frequency measurements appear to indicate many uncon- 
ventional properties of the Meissner and mixed states. 
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